Theory of indirect resonant inelastic X-ray scattering 
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We express the cross section for indirect resonant inelastic X-ray scattering in terms of an intrinsic 
dynamic correlation function of the system that is studied with this technique. The cross section 
is a linear combination of the charge response function and the dynamic longitudinal spin density 
correlation function. This result is asymptotically exact for both strong and weak local core-hole 
potentials. We show that one can change the relative charge and spin contribution to the inelastic 
spectral weight by varying the incident photon energy. 



Introduction. The probability for X-rays to be scat- 
tered from a solid state system can be enhanced by orders 
of magnitude when the energy of the incoming photons 
is in the vicinity of an electronic eigenmode of the sys- 
tem. Such resonant X-ray scattering experiments Q] are 
performed on e.g. the K-edges of transition metal ions, 
where the frequency of the X-rays is tuned to match the 
energy of an atomic Is-Ap transition, which is around 5- 
10 keV S H H H Q, @- At this resonant energy a 
Is electron from the inner atomic core is excited into an 
empty 4p state, see Fig. 1. In transition metal systems 
the empty 4p states are far (10-20 eV) above the Fermi- 
level, so that the X-rays do not cause direct transitions of 
the Is electron into the lowest 3<i-like conduction bands 
of the system. Still this technique is sensitive to exci- 
tations of electrons near the Fermi-level. The Coulomb 
potential of the Is core-hole causes e.g. very low en- 
ergy electron- hole excitations in the valence/conduction 
band: the core-hole potential is screened by the valence 
electrons. When the excited 4p-electron recombines with 
the Is core-hole and the outgoing photon is emitted, the 
system can therefore be left behind in an excited final 
state. Experimentally the momentum q and energy u> 
of the elementary excitation is determined from the dif- 
ference in energy and momentum between incoming and 
outgoing photons. Since the excitations are caused by 
the core-hole, we refer to this scattering mechanism as 
indirect resonant inelastic X-ray scattering (RIXS). 

RIXS is rapidly developing due to the recent increase 
in brilliance of the new generation synchrotron X-ray 
sources. High flux photon beams with energies that are 
tunable to resonant edges are now becoming available. 
At present an energy resolution of about 300 meV can 
be reached. In the near feature it seems experimen- 
tally feasible for RIXS to become sensitive to the low 
energy excitations of the solid, where excitation energies 
are of the order of room temperature. Such low- lying 
electronic excitations can, for example, be collective fea- 
tures such as plasmons, magnons, orbitons, excitons and 
single-particle-like continua related to the band struc- 
ture. RIXS provides a new tool to study these elementary 
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FIG. 1: Schematic representation of the indirect resonant in- 
elastic X-ray scattering (RIXS) process. 



excitations. 

For the interpretation of spectroscopic data, it is very 
important to express the scattering cross section for a 
technique in terms of physical correlation functions. It is 
well known, for example, that in standard optical experi- 
ments the charge density response function (or dielectric 
function) at q ~ is measured. In electron energy loss 
experiments this linear response function is measured at 
finite transferred momentum. 

In this paper, we determine the dynamical correlation 
function that is measured in indirect resonant inelastic 
X-ray scattering. For a local core-hole potential, the dy- 
namical correlation function turns out to be a linear com- 
bination of the charge density and longitudinal spin den- 
sity response function. We show that the actual linear 
combination that is measured depends on the energy of 
the incoming photons and we determine the precise en- 
ergy dependence of its coefficients. The expressions for 
the scattering cross section that we derive are asymptot- 
ically exact in the two limits of either a very strong or 
a very weak local core-hole potential. Our results, more- 
over, constitute a smooth interpolation between these 
two extremes. 

Series expansion of the scattering cross section. The 
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Kramers-Heisenberg formula [ft 1 1(31 for the resonant X- 
ray scattering cross section is 1 1 1 1 
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where / and i denote the final and initial state of the 
system, respectively. The sum is over all final states. The 
momentum and energy of the incoming/outgoing photons 
is q in /out and w? n/out and the loss energy u = w° ut - wf n 
is equal to the energy difference between the final and 
initial state ujfi = Ef — Ei. In the following we will take 
the energy of the initial state as reference energy: Ei = 0. 
The scattering amplitude Aft is given by 



A fi = 



(f\ 6\n)(n\6\i) 
- E n - iT ' 



(2) 



where Lo rcs the resonant energy, n denotes the interme- 
diate states and O the (dimensionless) dipole operator 
that describes the excitation from initial to intermedi- 
ate state and the de-excitation from intermediate to final 
state. The energy of the incoming X-rays with respect 
to the resonant energy is uj- ln (this energy can thus ei- 
ther be negative or positive: w- m = ujf n — oj tos ) and E n is 
the energy of intermediate state \n) with respect to the 
resonance energy. 

In the intermediate state a core-hole and a photo- 
excited electron are present. When we take the Coulomb 
interaction between the intermediate state core-hole and 
the valence band electrons into account, we obtain a fi- 
nite inelastic scattering amplitude. In that case there is 
a non-zero probability that an electron-hole excitation is 
present in the final state, see Fig.l. 

The intermediate state, however, is not a steady state. 
The highly energetic Is core-hole quickly decays e.g. 
via Auger processes and the core-hole life-time is very 
short. The Heisenberg time-energy uncertainty relation- 
ship then implies that the core-hole energy has an ap- 
preciable uncertainty. This uncertainty appears in the 
formalism above as the core-hole energy broadening T 
which is proportional to the inverse core-hole life-time. 
Note that the life-time broadening only appears in the 
intermediate states and not in the final or initial states 
as these both have very long life times. This implies that 
the core-hole broadening does not present an intrinsic 
limit to the experimental resolution of RIXS : the loss 
energy ui is completely determined by kinematics. 

When the incoming energy of the X-rays is equal to a 
resonant energy of the system u>i n ~ E n = and we see 
from Eqs. 11121) that the resonant enhancement of the X- 
ray scattering cross section is jofres/r) 2 , which is ~ 10 6 
for a transition metal K-edge |l2j. 

In a resonant scattering process, the measured sys- 
tem is generally strongly perturbed. Formally this is 
clear from the Kramers-Heisenberg formula Q , in which 



both the energy and the wavefunction of the intermediate 
state -where a potentially strongly perturbing core-hole 
is present- appear. This is in contrast with canonical op- 
tical/electron energy loss experiments, where the probing 
photon/electron presents a weak perturbation to the sys- 
tem that is to be measured. 

To calculate RIXS amplitudes, one possibility is to nu- 
merically evaluate the Kramers-Heisenberg expression. 
To do so, all initial, intermediate and final state energies 
and wavefunctions need to be known exactly, so that in 
practice a direct evaluation is only possible for systems 
that, for example, consist of a small cluster of atoms [l^ . 
In this paper, however, we show that under the appro- 
priate conditions we can integrate out the intermediate 
states from the Kramers-Heisenberg expression. After 
doing so, we can directly relate RIXS amplitudes to linear 
charge and spin response functions of the unperturbed 
system. For non-resonant scattering, one is familiar with 
the situation that the scattering intensity is proportional 
to a linear response function, but for a resonant scatter- 
ing experiment this is a quite unexpected result. 

Let us proceed by formally expanding the scattering 
amplitude in a power series 



A fi = 



iT 
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where we introduced the matrix elements 



Mi = 
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We denote the denominator of the expansion parameter 
dJi n — iT by the complex number A, so that 



= ^(f\d(H int ) l d\i), 



(5) 



where _ff; nt is the Hamiltonian in the intermediate state. 
We thus obtain the following series expansion for the res- 
onant cross section 
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d(w-w fi ). (6) 



Spinless fermions. We first calculate the resonant X- 
ray cross section in the case where the valence and con- 
duction electrons are effectively described by a single 
band of spinless fermions: spin, and orbital degrees of 
freedom of the valence electron system are suppressed. 
Physically this situation can be realized in a fully satu- 
rated ferromagnet. 

The final and initial states of the system are deter- 
mined by a Hamiltonian Hq that describes the electrons 
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around the Fermi-level. The generic form of the full 
many-body Hamiltonian is 



■3i 



(7) 



where i and j denote lattice sites with lattice vectors 
Hi and Rj and i,j range from 1 to N, where N is the 
number of sites in the system. The hopping amplitudes 
of the valence electrons are denoted by Uj and the c/c*- 
operators annihilate/create such electrons. The Coulomb 



interaction between valence electrons is Va 



|Ri-Rj 



as the Coulomb interaction only depends on the relative 
distance between two particles. 

The intermediate states are eigenstates of the Hamil- 
tonian H mt = Hq + H core , where H core accounts for the 
Coulomb coupling between the intermediate state core- 
hole and the valence electrons: 



lJ C j C 01 



(8) 



where Si creates a core- hole on site i. We assume that 
the core-hole is fully localized and has no dispersion. We 
will see shortly that this leads to major simplifications in 
the theoretical treatment of indirect RIXS. The core-hole 
- valence electron interaction is attractive: V c < 0. The 
dipole operators are given by 



= ]Te- iqi "- R < 



SiPi 



. e iqout-R is t p . 



+ h.c, 



(9) 



where p^ creates a photo-excited electron in a Ap state 
and h.c. denotes the Hermitian conjugate of both terms. 

In order to calculate the cross section, we need to eval- 
uate the operator {H- mt ) 1 in equation ©. A direct eval- 
uation of this operator is complicated by the fact that 
the initial and intermediate state Hamiltonians do not 
commute. We therefore proceed by expanding (-ffmt)' in 
a series that contains the leading terms to the scattering 
cross section for both strong and weak core-hole poten- 
tials. After that we will do a full re-summation of that 
series. We expand 



(H int yo\i) = (H + H core yd\i) 
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where we used [H ,O] — and H \i) = 0. For strong 
core-hole potentials the expansion is dominated by the 
m = term {H core f . Also for weak core-hole potentials 
we obtain the correct leading terms. It is straightforward 
to shown that for a strong core-hole potential the first 
correction to the expansion above is smaller by a factor 
t/V c and for a weak core-hole potential by a factor V c /t. 
Thus the expansion is well behaved as in the appropriate 
limits the leading order corrections are vanishing. 



Let us for the moment consider the strong core-hole 
potential limit and keep in the expansion only the term 
m = 0. Due to the fact that in the intermediate state 
only one core-hole is present and that this core-hole has 
no dispersion, we have 



M l (V c »t) = ±(f\Y / e^ R ' ^ 



vrAcA N>, (ii) 



where the transfered momentum q = q out — qi n • 

The first important observation is that the term I = 
does not contribute to the inelastic X-ray scattering in- 
tensity because M = (f\ J2i e 4q ' Ri |i) = Ndq^Sf^, which 
only contributes to the elastic scattering intensity at 
q = and other multiples of the reciprocal lattice vec- 
tors. From inspection of equation (0J we see immedi- 
ately that the I = term actually vanishes irrespective 
of the strength of the core-hole potential. This is of rel- 
evance when we consider the scattering cross section in 
the so-called " fast-collision approximation" [1J| . This ap- 
proximation corresponds to the limit where the core-hole 
life time broadening is the largest energy scale in system 
(r — > oo or, equivalently, im[A] — > — oo). In this limit 
only the / = term contributes to the indirect RIXS am- 
plitude and the resonant inelastic signal vanishes. In any 
theoretical treatment of indirect resonant scattering one 
therefore needs to go beyond the fast-collision approxi- 
mation. 

The second observation is that Mi is a 2'-particle cor- 
relation function. If we measure far away from reso- 
nance, where |i?e[A]| ^ 0, the scattering cross section 
is dominated by the 1 = 1, two-particle, response func- 
tion. When the incoming photon energy approaches the 
resonance, gradually the four, six, eight etc. particle re- 
sponse functions add more and more spectral weight to 
the inelastic scattering amplitude. Generally these multi- 
particle response functions interfere. We will show, how- 
ever, that in the local core-hole approximation the multi- 
particle correlation functions in expansion (|10|l collapse 
onto the dynamic two-particle (charge-charge) and four- 
particle (spin-spin) correlation function. 

In hard X-ray electron spectroscopies one often makes 
the approximation that the core-hole potential is local. 
This corresponds to the widely used Anderson impurity 
approximation in the theoretical analysis of e.g. X- ray 
absorption and photo-emission, introduced in Refs. |l5l 
Il6lll7fl . This approximation is reasonable as the Coulomb 
potential is certainly largest on the atom where the core- 
hole is located. 

We insert a local core-hole potential V£ = USij and 
find after performing the sum over m in equation (|10fl 
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Using that £V e lq ' Ri c- c 4 = X)k c k+ q c k = Pq is tne den- 
sity operator, we now perform the sum over / in equa- 
tion and obtain 



A% = Pi(uj,U)(f\ P(i \i) 



(13) 



where the resonant enhancement factor is P\{oj,U) = 
Uui les [(A — U)(A — u})]^ 1 and u — Ej. For spinless 
fermions with a local core potential the indirect RIXS 
cross section thus turns out to be the density response 
function -a two-particle correlation function- with a res- 
onant prefactor Pi(w) that depends on the loss energy u>, 
the resonant energy lu tcs , on the distance from resonance 
cji n (= i?e[A]), on the core-hole potential U and on the 
core-hole life time broadening T(— —Im [A]). The den- 
sity response function is related to the dielectric function 
e(q, u>) and the dynamic structure factor S(c[, w) [18|, so 
that for fixed U we obtain 



dildtj 



oc -|PiH| 2 Im 
ex |PiH| 2 S(q,uj). 



1 



U q e(q,w) 



(14) 



Fermions with spin. We generalize the calculation 
above to the situation where the electrons have an addi- 
tional spin degree of freedom. In the Hamiltonians (|7I8(1 
we now include a spin index a (with a =| or J.) to the an- 
nihilation and creation operators: c; — > Ci a and Cj — ► Cj>/ 
and sum over these indices, taking into account that the 
hopping part of the Hamiltonian is diagonal in the spin 
variables. In order to re-sum the series in equation © 
we now need to evaluate expansions of the number oper- 
ators of the kind (n-\ + n^) 1 . For / > 0, this leads to 

A fi = [P 1 H-P 2 H](/|S q |z)+P 2 H(/| Pq | l ), (15) 

withP 2 (w,C7) =Pi(2w, 2C/)/2 andp q = p q +p q . Here we 
also introduced the longitudinal spin density correlation 
function S q = igrg+fi J2k ^k+ q ■ S_k- Clearly the contri- 
butions to the scattering rate from the dynamic longitu- 
dinal spin correlation function and the density correlation 
function need to be treated on equal footing as they inter- 
fere 01 . Moreover, the spin and charge correlation func- 
tions have different resonant enhancements. For instance 
when Pe[A] = U, the scattering amplitude is dominated 
by Pi(u>) and hence by the longitudinal spin response 
function. At incident energies where Pe[A] = 2U, on the 
other hand, P 2 (u>) is resonating so that the contributions 
to the inelastic scattering amplitude of charge and spin 
are approximately equal. 

Note that we only considered single band systems in 
this Letter. As a consequence the strength of the core- 
hole Coulomb potential is the same for all electrons in 
the valence band. In multi-band systems |2(ij . however, 
the Coulomb coupling will depend on the actual orbitals 
and bands that the electrons are in. For electrons in 
localized, atomic-like orbitals (e.g. transition metal 3d 



orbitals) the Coulomb interaction with the core-hole is 
largest. Indirect RIXS is therefore mostly sensitive to 
excitations of these localized electrons. 

Conclusions. We presented a series expansion of the 
indirect resonant inelastic X-ray scattering amplitude, 
which is asymptotically exact for both small and large 
local core-hole potentials. By re-summing the series, we 
find the dynamical charge and spin correlation functions 
that are measured in RIXS. The calculation of these cor- 
relation functions from e.g. model Hamiltonians or by 
means of ab initio electronic structure methods, will allow 
for a direct comparison of computed response functions 
with experimental RIXS spectra. 
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